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Abstract 


Two  theoretical  investigations  are  given  relating  to  SAR  imaging  of 
perturbed  surface  waves.  The  first  is  an  analysis  of  the  effect  of  a  long 
wavelength  corrugation  on  the  SAR  image  formed  by  lstrorder  Bragg  scattering 
from  a  deterministic  sea  wave  pattern  such  as  a  Kelvin  wake.  A  standard  SAR 
processor  is  assumed.  vThe  second  is  a  stochastic  treatment  of  the  SAR 
imagery  formed  by  a  lst^order  Bragg  scattering  process  from  surface  waves 
modulated  by  the  horizontal  surface  currents  of  internal  waves. 
Propagational  velocities  (normalized  by  the  SAR  platform  velocity)  are 
included  to  first  order.  For  the  Kelvin  wake  analysis  it  is  shown  that  the 
vertical  curvature  of  the  corrjigating  wave  can  rotate  the  effective  Bragg 
acceptance  angle  by  up  to  for  typical  wave  parameters.  For  the 
internal  wave  analysis,  SAR  resolutions  are  given,  and  in  particular  it  is 
shown  that  typical  azimuthal  velocity  components  are  important  in 
determining  the  azimuthal  resolution  dimension  of  the  scene,  especially  for 
resolutions  of  a  few  metres  size.  ^ 


RAsumA 


Le  texte  dAcrit  deux  Atudes  thAoriques  sur  1' imager ie  SAR  de  vagues 
de  surface  perturbAes.  La  premiAre  est  une  analyse  de  l'effet  d'une 
cannelure  de  grande  longueur  d'onde  sur  1' image  SAR  formAe  par  la  diffusion 
de  Bragg  du  premier  ordre  A  partir  d'un  diagramme  dAterministe  des  vagues  de 
1'ocAan,  par  exemple  un  sillage  Kelvin.  On  suppose  le  recours  A  un 
processeur  SAR  standard.  La  deuxiAme  Atude  consiste  en  un  traitement 
stochastique  de  l'imagerie  SAR  formAe  par  la  diffusion  de  Bragg  du  premier 
ordre  de  vagues  de  surface  modulAes  par  les  courants  de  surface  horizontaux 
des  ondes  internes.  Les  vitesses  de  propagation  (normalisAes  par  la  vitesse 
de  la  plate-forme  SAR)  sont  donnAes  jusqu'au  premier  ordre.  L' analyse  de 
sillage  Kelvin  perraet  de  dAmontrer  que  la  courbure  verticale  de  la  vague  de 
cannelure  peut  faire  tourner  jusqu'A  10°  1' angle  d' acceptation  efficace  de 
Bragg  pour  des  paramAtres  de  vagues  typiques.  Dans  le  cas  de  1' analyse  des 
vagues  internes,  on  prAcise  les  valeurs  de  la  rAsolution  SAR,  et  on  constate 
plus  particuliArement  que  les  composantes  de  vitesse  azimutale  typiques  sont 
importantes  pour  dA terminer  la  dimension  de  rAsolution  azimutale  de  la 
scAne,  en  particulier  pour  les  rAsolutions  de  1' ordre  de  quelques  mAtres. 
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1 .  Introduction 


Several  studies1,2,3  have  been  made 'of  the  effect  of  motion  in  scene 
elements  on  SAR  processing.  In  those  studies  the  nature  of  the  motion  has 
been  either  random1,  in  which  case  temporal  coherence  estimates  have  been 
given  with  their  effect  on  the  final  image.,  or  organized  in  the  form  of 
surface  waves2  but  with  random  phasing  relative  to  the  Bragg  scatterers.  In 
either  case  it  has  been  shown  that  "velocity  bunching"  takes  place,  i.e., 
image  misregistrations  occur  arising  from  misinterpretations  of  phase 
histories  due  to  the  combination  of  platform  velocity  and  scene  element 
velocity.  The  decorrelation  time  associated  with  random  motions  also  serves 
to  increase  the  effective  azimuthal  SAR  resolution  distance3.  Scene 
acceleration  effects  can  also  be  important2  producing  direct  changes  in 
resolution. 

Two  separate  problems  are  discussed  in  the  present  paper,  SAR 
imagery  of  (a)  Kelvin  wakes,  and  (b)  moving  internal  wave  fields.  In  both, 
a  spectral  description  is  used  for  the  Bragg  scatterers  rather  than  a 
correlation  description.  This  offers  a  simplified  interpretation  based  on 
narrow-band  processes.  For  the  Kelvin  wakes,  a  deterministic  problem  is 
solved  in  which  the  Bragg-scattering  field  exists  (the  diverging  wake  wave- 
field)  in  the  presence  of  a  long  "corrugating"  wave  (the  transverse  wake 
wave-field).  The  Kelvin  wake  is  assumed  to  be  steady,  and  so  the  Bragg- 
scatterers  are  phase -locked  to  the  corrugating  wave:  specifically,  both  have 
the  same  component  of  phase  velocity  along  the  ship's  track,  namely,  the 
ship's  speed.  The  corrugating  field  is  not  particularized  to  a  detailed 
Kelvin  wake  structure  but  is  left  as  an  infinite  sinusoid.  The  Bragg 
scatterers  are  also  left  unspecified  in  detail.  The  problem  is  solved  in 
terms  of  a  filter  function  which,  when  applied  to  the  Bragg  scattering 
divergent  wave  field,  provides  the  unsquared  complex  SAR  image.  A  numerical 
example  is  given  for  conditions  pertinent  to  narrow-V  wakes  recorded  in 
Dabob  Bay ,  WN4 . 
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made,  and  for  the  present  purpose  of  discussing  lst-order  Bragg  backscatter, 
only  terms  linear  in  surface  wave  height  will  be  kept.  The  ocean  will  also 
be  treated  as  being  a  perfect  conductor  and  only  horizontal  polarization  in 
the  incident  radar  field  will  be  examined.  These  last  two  conditions  can  be 
removed  but  only  with  some  additional  algebraic  complexity  that  does  not 
illuminate  the  central  issues  of  the  present  problems.  The  source  and 
receiver  will  also  be  modelled  as  points,,  but  with  beam  shaping  introduced 
at  a  later  stage. 

With  these  approximations,  the  scattered  field  becomes 


$(x,z,t) 


z2  ff  .cR)3"(ts)  d£< 

4tt2c2  -oo  tl2r22 


(1) 


where  is  the  vector  radar  field,  z  is  the  source  -  receiver  altitude, 

x(=»x,y)  is  a  coordinate  frame  located  at  and  moving  with  the  source - 
receiver,  c  is  the  speed  of  radar  propagation  in  air,  is  represents  the 
radiated  radar  field  (with  units  of  ^’length),  and  x'(=x',y')  is  a  non¬ 
moving  coordinate  frame  fixed  to  the  ambient  ocean  surface.  Furthermore, 


tR  -  t  -  ri/c  ,  (2) 

ts  -  tR  -  r2/c  -5L-  (y'-vt)  “  t  -  £il£2  ,  (3) 

c2  c 

ri  -  J (x-x')2  +  (y+Vt-y')2  +  z2  ,  (4) 

r2  -  Jx'2  +  (y'-VtR)2  +  z2  -  Jx'2  +  (y'-Vt)2  +  z2  ,  (5) 

2"(ts)  -  d2S(ts)/dts2  ,  (6) 


T  -  surface  waveheight  , 
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and  here  V  is  the  aircraft  speed  oriented  along  the  y-direction  (azimuthal) , 
and  the  centre  of  the  radar  beam  is  oriented  along  the  x-direction  (range) . 
The  ratio  V/c  is  sufficiently  small  that  the  approximate  forms  of  Equations 
(3)  and  (5)  will  be  used. 

The  rx-  and  r2 -terms  in  the  denominator  of  the  integrand  are 
sufficiently  slowly  varying  and  sufficiently  similar  that  they  can  be 
removed  from  the  integral  and  each  set  equal  to  R,  the  distance  from  the 
source-receiver  to  the  centre  of  the  scene  element  (i.e.  pixel).  !>  is  also 
narrowband  enough  in  frequency  that  the  approximation 

2"  -  -uPt  (7) 

can  be  used  with  no  appreciable  loss  of  accuracy  for  the  present  purposes, 
with  wc  being  the  central  frequency  in  the  radiated  radar  field. 

The  vector  character  of  $  and  3  appears  only  in  the  radiated  fields 
and  it  is  scalar  versions  of  these  that  are  used  in  the  SAR  processing.  The 
conversions  take  place  at  the  antenna  and  bring  in  the  beam- shape  function. 
The  received  voltage  is  given  by 


CO 

•  v(t)  -  ///  fR(x,y,z).^(x,yIz,t)dxdydz 
-00 


-  Iff  FR(x,y,z) 
-00 


x  JT  ///  — — ~'Rl-S  (yV*s’yS,zs>  dx'dxsdysdzsdxdydz  (8) 
-oo  ri'i(x,y1z,t,x')r2^(xSlys,z0,t,x') 
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Here,  ?T  is  the  beam-shaping  function  for  transmission,  ?R  is  the  beam¬ 
shaping  function  for  the  receiver,  and  the  vector  character  of  the  radiated 
field  is  included  explicitly  in  ?T  with  Si  now  merely  being  the  (scalar) 
time -dependent  part.  Also,  the  receiver  altitude  is  now  at  z  rather  than 
being  coincident  with  the  transmitter  at  zs.  Equation  (8)  can  be  rewritten 
as 


H t)  -  fff  ?R(x,y,z)  ///  ?T(xs,ys,zs)<£p(x,y,z;xs,ys,z0;t)dxdydzdxsdysdzs  (9) 


with 


( t) 


22s  rf  g(x‘ ,tR)Sj(ts)  3X, 
4tt2c2  .a,  ri2r22 


(10) 


Because  of  the  far- field  condition  for  the  scattering  and  all  the 
preceding  assumptions ,  Equation  (9),  with  Equation  (10)  can  be  simplified 
into 


0(t)  =  - 


z2ca^ 

47t2c2rA 


If  t(x’ 


tRjSCtso''*'  .y'  *vt,z)dx' 


(11) 


where 


S(t;x,y,z)  -  S1(t)B(x,y , z)  (12) 

and  B  is  now  the  modulus  squared  of  the  voltage  beam- shape  factor.  In  this 
form,  beam  shape  is  incorporated  directly  into  the  scalar  S  which  describes 
the  incident  radiation. 

The  SAR  image  is  formed  by  matched- filtering  ^(t)  with  a  version  of 
4>  that  would  be  obtained  from  a  stationary  point  reflector.  The  matching  is 
done  over  the  range  coordinate  i.e.,  over  a  t- interval  that  corresponds  to 
the  radar  field  propagation  time  from  the  antenna  to  the  most  distant 
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scatterer  to  be  imaged  and  back  to  the  antenna.  Matching  is  also  done 
separately  over  azimuth  and  for  that  S  and  4>  are  modelled  as  being  pulse  - 
trains  and  the  matching  summation  is  performed  over  a  suitable  number  of 
contiguous  pulses,  a  number  which  can  represent  an  aircraft  travel  distance 
that  is  a  significant  fraction  of  the  azimuthal  beamwidth  at  the  sea 
surface . 


For  an  explicit  pulse- train  description,  Equation  (11)  becomes 


Mr) 


zi?c2  //  9 (x ,  tn+r -  £iS)S(r  -  IllljEln  ;  x' ,y' -Vtn, zQ)dx' 
4*-2c2r4  _ro  c  c 


(13) 


where  r  is  time  measured  from  the  start  of  the  nth  pulse,  t^,  and  rln  and  r2n 
are  ranges  at  t„.  For  most  applications  r  can  be  dropped  from  the  time- 
argument  of  9  because  the  radar  scattering  interval  is  effectively  only  a 
small  fraction  of  a  millisecond,  a  time  short  enough  that  the  surface  height 
usually  does  not  change  appreciably.  The  same  simplification  has  been  made 
in  the  spatial  argument  of  S  (relating  to  the  beam  shape  factor  in  azimuth) , 
because  the  aircraft  travels  a  very  short  distance  in  the  radar  propagation 
time  and  the  beam  shape  is  smooth  enough  that  this  distance  is  unimportant. 


The  point-reflector  replica  of  <£(t)  is  obtained  by  letting  9  be  a 
fixed  5-function  at  location  (x0,y0).  The  coordinate  y0  is  left  unspecified 
for  the  present  but  it  will  eventually  be  set  to  zero  because  the  processing 
will  always  be  considered  to  be  done  symmetrically  in  azimuth  about  the 
location  of  the  replica.  At  this  stage  also  the  aircraft  speed  mismatch  can 
be  included  by  using  Vx  rather  than  V  in  the  replica  description,  giving 


^rep ( T ) 


2^c2°o  S(r  .  2jx024-(y0-V1tn)2+z02, 


4^2c2d  4 


(14) 


Here,  cr0  is  the  scattering  strength  of  the  point  reflector. 
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For  a  Gaussian  beam  shape  and  Gaussian  pulse  shape,  S  takes  the 


form 


S(t;x,y-Vtn,z)  =  e 


-io)ct-iw't2-t2/TG-  (y-Vtn)2/Yc 


+  c  .  c  . 


(15) 


where  TG  and  YG  are  the  Gaussian  bearawidth  factors  in  time  and  azimuth 
respectively  and  where  c.c.  indicates  complex  conjugate.  Here  also, 
frequency  modulation  is  explicitly  included  by  means  of  the  w'-term,  and  the 
range  coordinate  beam- shape  is  ignored  because  Lt  is  usually  very  broad  and 
contributes  effectively  only  to  a  general  shading  in  the  overall  image. 


In  carrying  out  the  matching  process,  <f>n(r)  and  <£rep(r )  are  both 
demodulated  ( frequency- shifted  to  base -band)  and  only  one  side  of  the 
transform  is  used  (e.g.  positive  frequencies  only).  This  results  in  some 
computational  simplicity  for  digitally  formed  images  although  neither  step 
is  crucial  to  image  formation.  Appropriate  range -dependent  phase  -  shifting 
terms  are  also  incorporated  to  eliminate  extraneous  modulation  products  in 
the  result.  The  undetected  image  is  given  by 


W(x0)  -  Z  J  *n('>**  (Odr  (16) 

n  reP 

-  00 

where  the  <f>' s  are  now  modified  as  just  described  and  the  sum  is  over  N 
contiguous  pulses. 

By  using  Equations  (13)  and  (14),  Equation  (16)  can  be  written 


rl',  2 


W(xQ)  =  S  ff  ?  (x'  ,  tn  -  I±I!)  f  (  -  Z-^P-  So(r  -  rln+r2n;x>  >y-  -Vtn,z) 

C  .00  47r2C2RA  c 


n 


2io)cx02 


_  2lxV.^yi'=n.)2ii,.2;x.iy..VltniZo),ecJ^02 


dx’dy'  (17) 
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where  the  subscript  D  on  each  S  indicates  demodulation,  and  positive 
frequencies  only,  and  the  unimportant  amplitude  scale  factor  term  in  $rep, 
i.e.  z02wc2<70/4jt2c2R04  ,  has  been  dropped.  The  purely  phase- shifting  term 
exp{2ia)cx0  */JZ  2+z02 )  is  also  shown  explicitly  (although  it  need  not  be 
inside  the  integral) . 

For  analytical  simplicity,  the  sum  on  n  is  approximated  as  an 
integral  over  t„  using  the  definition 

tn  -  nA  (18) 

where  A  is  the  time  interval  between  successive  radiated  pulses.  With  this 
transformation,  the  sum  becomes 

2  -  J  /  ...  dtn  .  (19) 

One  further  simplifying  change  can  be  made,  shifting  from  tn  to  the  time- 
argument,  t' ,  of  ?  ,  i.e. 

t'  -  <Vrln/c  .  (20) 

To  the  desired  degree  of  accuracy,  and  using  an  impulse  response 
formulation,  Equation  (17)  becomes 

W(x0)  -  J/J  t(x' ,t')I(x' ,x0,t')dx'dt' 


-CO 


(21) 
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where 


I(x,x0,t' ) 


-z2g,c2 

4w2c2r^a 


exp 


2iav 


,  x02+z02j 


►  J  Si)(r  - 


£I±£2;  x'.y'-Vt' 


X  S *(r  -  2Jxo2+<yolYl-t.'  -Vlr.l/.cl2tl°? ;  x'.y'-Vit'-  XlEl|ZO)dr 
Uc  e 


(22) 


and 

-  iw'  (t-e)2  -  &:02  -  Yi_ 

SD(t-« ;x,ylz)  -  e  c  TG2  YG2  (23) 

The  r's  in  Equations  (20)  and  (22)  are  actually  functions  of  t  (see 
Equations  (4)  and  (5))  but  to  the  required  accuracy  they  may  be  treated  as 
time- independent ,  particularly  as  used  in  Equation  (22).  This  can  be  shown 
by  expanding  tR  and  t  about  t,,  and  substituting  these  expansions  into  the 
expressions  for  rx  and  r2  and  ignoring  terms  of  order  V/c  that  are  small. 

With  all  the  approximations  incorporated,  the  final  expression  for 

I  becomes 


I(x,x0,t')  -  - 


z^ 


4w2c2r4a 


exp 


2iwr 


+i|c  (r1+r2-2r0) 


+Zr 


/(y'-Vit')2+(y'-Vt')2 

l  .  (r1+r2-2r0)2  2  * 

^  YG2 

1  2c2tg2  G  J 

(24) 


where 


rl 


2+(y'-Vt') 


2+z2 


(25) 


r2  “  rl 


(26) 
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and 

r0  “  Jx02+(y0-Vit' )2+z02  .  (27) 

This  completes  the  general  analysis.  Fourier  transform 
descriptions  for  <r  will  be  used  in  the  following  sections  but  with  forms 
that  are  specifically  relevant  for  each  problem.  Also,  the  r's  of  Equations 
(25)  to  (27)  will  be  binomially  expanded  about  a  coirmon  position  and  terms 
up  to  and  including  second  order  in  deviations  from  this  position  will  be 
kept.  For  the  Kelvin  wake  problem,  z  in  Equation  (25)  will  also  be  allowed 
to  be  a  function  of  position  in  order  to  describe  the  effect  of  the  height 
of  the  corrugating  wave . 

3 .  Kelvin  Wakes 


The  entire  wake  wave -field  including  the  corrugating  portion  is 
reducible  to  a  steady- state  by  a  simple  transformation  to  coordinate  axes 
moving  vith  the  source  ship.  Thus,  with  cs  as  the  ship's  velocity,  and  with 
distance  measured  from  the  ship  being  given  by  rj ,  the  following  pertain: 

V  ~  x’ -cst'  ,  (28) 

C(x'.t')  -  c07)  , 

and,  from  Equation  (21) 

W(xQ)  “  //  ?(»7 )  /  I(»f+cst'  ,x0,t'  )dt'd?7  (29) 


The  divergent  wave- field  c  can  now  be  expressed  as  a  two-dimensional  Fourier 
transform, 

«(*>  "  J7  P(ic)ei^r7dij 


-00 


(30) 
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so  that 


p(lc)  -  //  dr?  , 

4jt  _  oo 


(31) 


and  W(xQ)  -  JJ  p(Ic)e^*xo  (J/J  eik*(r7-x0)  I (r?+cst '  , xQ  ,  t '  )dt '  dr^ ) die  (32) 


The  term  eik*xo  taken  out  explicitly  in  the  integrand  of  this  last 

expression  in  order  that  W(x0)  may  be  seen  to  be  simply  a  filtered  version 
of  ?  by  comparing  Equations  (32)  and  (30).  The  filter  function,  denoted  Q2, 
is  given  by  the  term  in  the  braces  in  Equation  (32),  i.e., 


Q2(Ic;x0)  -  /JJ  ei^(»?-xo)l(^+cst',x0,t')dt'd^  (33) 

-  CO 

and  it  contains  all  of  the  frequency- shaping ,  phase  -  shifting  and 
demodulation  that  the  SAR  process  applies  to  the  <;  -field  in  the  Fourier 
domain.  In  this  steady-state  formulation,  all  of  the  fluid- dynamically 
interactive  effects  of  the  corrugating  field  and  its  currents  on  the  ?- field 
are  contained  in  p(ic)  and  need  not  be  expressed  explicitly.  This  simple 
situation  will  not  hold  in  the  next  section  which  deals  with  internal  wave 
modulations  that  have  a  stochastic  nature  and  there  the  frequency- shifting 
effect  of  the  currents  in  the  corrugating  field  will  need  to  be  explicitly 
incorporated. 

To  progress  further  with  the  Kelvin  wake  problem  it  is  necessary  to 
return  to  Equations  (24)  to  (27)  and  expand  the  expressions  for  rltr2  and  r0 
using  the  definitions 


x'  -  x0  +  M  . 

(34) 

y'  -  Vt'  +y"  , 

(35) 

y0  -  Vit'  +  q"  , 

(36) 

z  -  zQ  +  tT(rT)  , 

(37) 

Rq  ”  Jx02+z02 

(38) 
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With  terms  expressed  up  to  and  including  second  order  in  m.  y",  q" 


and  ?T, 


R  +  Mxo+  ?tzq  +  u2+y"2+9T2 


xo242  .  M?TXqZq 


-  zo2^T2  + _ 


r°  “  R°  +  2R^  +‘ ‘ • 


and,  from  Equations  (24)  and  (34)  to  (36) 


Kx'.xo.t')  -r  I(M,y" ,x0,q") 


zo^c4  « 
4tt2wc2R0^A 


r  2iojc  Xo2  +  „  qll)  .  y"2+(y"  +  (V-Vj) Al(yQ-q")  )2 

c  R0  c  ’  ’  Yg2 


-  (/ix0+?Tz0) 2  («'  2TG^- 1) /2c2TG2R02 


where 


*(A‘.y".q")  -  rx+r2-2rc 


-  ?AiXo+?Tzo  +  "2zo2  +  v''2-an2  +  ?T2xo2  .  2^?Txozo 


and  where,  by  Equations  (28),  (34)  and  (36) 


A*  -  »?x-xo+Csx(yo-q")/vl  , 
y"  -  ^y+(csy-V)(y0-q")/V1 
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Here,  tjx  and  riy  are  the  components  of  r) ,  and  cax  and  csy  are  components  of 
-► 

cs- 

The  integration  in  t'  in  Equation  (33)  can  be  performed  after 
converting  t'  to  qM  with  the  use  of  Equation  (36),  and  the  result  is 


where 


Q2<£;x0)  -  J7  eik*  (»7 -xo)K(»7  ,x0)d»? 
-00 


(46) 


K(^,x0) 


zo2(Jc2  J[TG_eAl+Bi2/4Ci 
47r2c2R04AV1  JYjciT 


(47) 


To  obtain  this  expression,  the  argument  of  the  exponential  in  Equation  (42) 
has  been  written  as 

arg  -  Ai  +  Biq"-Ciq"2  (48) 

and,  by  substitution  of  Equations  (44)  and  (45)  into  (43)  and  (42)  and 
subsequent  inspection,  the  coefficients  Ax ,  Bx  and  Cx  can  be  readily 
obtained.  They  will  not  be  expressed  explicitly  here,  but,  because  of  the 
original  quadratic -form  expansions  of  rx,  r2  and  r0,  it  can  be  noted  that 
arg  itself  is  at  most  quadratic  in  products  of  rf  and  q" ,  and  so,  Ax  is 
quadratic  in  rj y  Bx  is  linear  in  rj  and  Cx  is  constant.  The  same  orders  exist 
for  sT.  A  further  simplification  will  also  be  made  to  these  coefficients, 
namely  that  terms  quadratic  in  the  small  quantities  |cs|/V  and  ?T/R0  will  be 
ignored.  For  the  usual  range  of  SAR  parameters,  these  quantities  are  of 
order  0.1  or  less  and  the  ignored  quadratic  terms  are  always  dominated  by 
linear  counterparts. 


Without  specifying  the  functional  form  of  ?T  on  rj ,  the  Fourier 
transform  represented  by  Equation  (46)  cannot  be  carried  out.  However, 
insight  can  be  gained  into  possible  effects  of  ?T  in  the  resulting  image,  by 
examining  the  special  case  in  which 


(49) 
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i.e.,  the  source  ship  travelling  parallel  to  the  SAR  aircraft  track,  and 

?T  “  ?o  +  +  l'72y(^r)o  •  (50) 

(For  simplicity,  the  processor  mismatch  is  also  set  equal  to  zero,  i.e., 
V-Vx.) 


The  surface  underlying  the  Bragg-wave  is  thus  given  a  tilt  and  a 
curvature  in  the  azimuthal  direction  only.  The  argument  of  the  exponent  in 
Equation  (47)  can  be  put  in  a  quadratic  normal  form 


arg  -  A  -  [  (*7X-Xs)/Resx]2lnl6  -  [  (r?y-ys)/Resy  ]  2lnl6  (51) 

where  Resx  and  Resy  represent  full -widths  at  half -maximum  and  xs  and  ys 
represent  the  position  in  tj  at  which  the  maximum  occurs.  For  the  case 
presented  by  Equations  (49)  and  (50), 


Resx 


?R0  Jl^4 

w'XqTg _ 

{1+1/w'2tg^} 1/2 


(52) 


Jln4 

Resy  -  - 2sIS - - - (1  +  £SX)  (53) 

(l-4(«^fi)2(*a)2Zo(«£4)  ,1/2 

wc  yG  5y2  ° 

The  expression  for  Resx  is  unmodified  by  the  presence  of  <;T  (to  this  order)  , 
but  it  can  be  seen  that  Resy  can  be  significantly  altered  by  cT's  curvature. 
This  gives  rise  to  the  possibility  of  brightened  spots  appearing  in  a  wake 
image  at  positions  where  there  is  a  coincidence  of  maximum  (positive) 
a  ti/dy  i  and  Bragg  scatterers  in  the  diverging  wake  components.  For 
stronger  conclusions  the  form  of  ?T  will  be  specified  more  fully  and 
numerical  evaluations  will  be  performed. 
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The  specific  values  chosen  for  numerical  treatment  pertain  to  the 
Canada  Centre  for  Remote  Sensing  L-band  system  (originally  constructed  by 
the  Environmental  Research  Institute  of  Michigan)  as  configured  in  July  1983 
and  used  in  overflights4  of  Dabob  Bay,  WN.  For  JOWIP  8/2 


uc  -  2w  (1 . 185) 109  rad/s  (54a) 

w'  -  2n  (33 . 3) 1012  rad/s  (54b) 

xQ  -  3033.5  m  (54c) 

z0  -  7010.0  m  (54d) 

yQ  -  0  m  ( 54e) 

Yg  -  162.1  m  ( 54f ) 

Tg  -  0.951  10'6  s  ( 54g) 

where  the  last  two  values  are  Gaussian  equivalents  which  result  in 

unperturbed  resolution  widths  of  ,j"20  m  in  both  directions.  For  JOWIP  8/2, 

V  -  134  m/s  (54h) 

cs  -  8  m/s  along  a  direction  -13°  off  azimuth  (54j) 


and  the  amplitude  of  the  corrugating  wave  had  a  value  of 

?0  “  0.1  m  (54k) 

at  a  distance  of  -2.5  km  behind9  the  source  ship  (USS  QUAPAW) .  For  this 
ship  speed,  the  wavelength  of  the  corrugating  transverse  wake  component  is 

As  -  41  m  (54m) 


calculated  from  the  gravity  wave  formula  Inc^/g. 
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The  corrugating  wavefield  ?T  is  modelled  as  a  simple  sinusoid. 
This  idealization  is  used  because  variations  in  actual  cT's  perpendicular  to 
the  wake  axis  are  minimal,  particularly  near  the  axis  where  the  L-band  Bragg 
waves  are  located  in  the  diverging  field.  Thus, 


<TT  “  sDcos  [2£(f7xsin0+r?ycos5  )+<f>)  (55) 
^  s 

where  9  is  the  angle  of  the  ship's  track  with  respect  to  the  aircraft  track 
(13°  for  the  present  example)  ,  and  4>  is  a  parameter  that  allows  the  full 
range  of  phase  of  ?T  to  be  investigated  numerically. 

In  the  Fourier  definition  of  Q2,  Equation  (46),  the  ^-coordinate 
system  can  be  rotated  through  the  angle  9  so  that  ?T  becomes  a  function  of 
only  one  variable.  Integration  can  be  performed  along  the  orthogonal 
component  of  the  rotated  rj  system,  and  since  the  exponential  argument  of  K 
is  a  quadratic  form  in  rj ,  apart  from  cT,  it  is  also  a  quadratic  form  in  the 
rotated  system,  and  the  integration  simply  returns  a  counterpart  quadratic 
form.  The  remaining  integral,  along  the  direction  9  must  be  done 
numerically. 

In  the  absence  of  a  corrugating  perturbation,  Q2  is  a  general  two- 
dimensional  Gaussian  form  in  £  and  it  can  be  described  by  its  maximum 
height,  the  wavenumber  position  at  maximum,  and  the  widths  along  the  U-axes 
to  the  half-maximum  points.  This  is  illustrated  in  contour  form  in 
Figure  1.  If  the  peak  is  also  rotated,  the  angle  of  rotation  is  necessary 
to  complete  the  description.  In  the  presence  of  the  perturbation,  these 
parameters  change,  and  the  effect  of  the  perturbation  can  be  described  fully 
by  the  amounts  of  the  changes.  Figure  2  shows  the  changes  in  form  of  the 
half -maximum  contour  of  Q2  for  different  phases,  i.e.  different  values  of  <f> 
in  Equation  (55),  along  the  corrugating  wave.  Figures  1  and  2  both  use  the 
parameters  specified  in  Equations  (54a)  to  (54m)  except  for  Equation  (54k). 
Figure  1  uses  <ro-0  and  Figure  2  uses  ?o-0.5  m. 
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Figure  1.  Contours  of  Q2  in  the  absence  of  a  corrugating  wave  (i.e.  ■?<,-()) 
but  with  the  source  travelling  along  5-13°.  The  lst-Bragg 
wavenumber  for  this  case  is  ( -19 . 712 -0 . 13238  .. )  rad/m.  The 
dotted  line  indicates  the  extreme  in  propagation  direction 
(2.20°)  for  Q2  at  its  half -maximum. 

KY  (rad/m) 


•5.0  -4.0  -i0  -20  -10  0.0  to  10  10  *0  5.0 


Figure  2.  The  half-maximum  contour  of  Q2  for  various  phases  (degrees)  along 
a  corrugating  wave  of  amplitude  ?o-0.5  m  propagating  at  5-13°  to 
the  azimuthal  direction.  Note  different  scales  from  Figure  1. 
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Without  a  corrugating  wave,  it  can  be  seen  from  Figure  1  that 
surface  waves  with  a  wavenumber  centered  at  - 19 . 712 . . , -0 . 132  rad/m  and 
within  a  circular  band  of  width  ±0 . 6240 ...  rad/m  are  accepted  by  the 
processor  with  a  gain  of  between  1  and  0.5.  This  bandwidth  implies  that  all 
waves  propagating  at  an  angle  within  ±1.81°  centered  at  0.38°  off  this  range 
direction  and  with  a  kx-value  of  -19.712  rad/m,  are  accepted  by  the 
processor  with  a  gain  reduction  of  no  more  than  a  factor  of  2.  (The  slight 
offset  of  the  peak  from  ky-0  is  due  to  the  non-zero  value  of  8,  i.e.  the 
non-zero  value  of  cS3t.  This  range  -  component  of  motion  in  the  Bragg 
scattering  wavefield  coupled  with  the  aircraft  motion  and  time  interval 
between  successive  radar  pulses  results  in  a  slight  skewing  of  the  Bragg- 
wave  field  in  the  SAR  pattern.) 

With  a  corrugating  wave  as  used  in  Figure  2,  it  can  be  seen  that 
the  surface  wave  wavenumbers  accepted  by  the  processor  vary  over  a  much  more 
considerable  range,  allowing  waves  propagating  at  up  to  11.5°  to  the  range 
direction  to  be  "seen"  (at  <£-90°)  ,  with  a  gain  reduction  of  no  more  than  a 
factor  of  2  from  the  peak  value .  This  variation  is  depicted  more  fully  in 
Figure  3  where  the  extreme  ky- value  for  the  half -maximum  contour  is  shown  as 
a  function  of  <j>  for  ?0-values  from  0  to  0.5  m.  In  a  similar  fashion, 
extreme  k^-values  are  shown  in  Figure  4,  and  the  Q2- value  at  its  peak  is 
shown  in  Figure  5,  for  parameters  as  in  Figure  3. 

From  Figure  3  it  can  be  seen  that  the  extreme  ky-values  are  largest 
for  <£“110°  (and  smallest  at  290°).  Using  the  central  k^-value  of  -19.712 
rad/m,  the  extreme  values  represent  surface  waves  propagating  at  angles  with 
respect  to  the  range  direction  as  shown  in  Figure  6.  Here  the  propagation 
angles  are  given  as  a  function  of  ?0,  and  propagation  directions  that 
pertain  to  the  center  of  the  processor  peak  are  also  given.  The  displayed 
data  shows  strongly  linear  dependence  with  ?0,  and  very  little  variation  for 
8 -values  over  the  range  9-0°  to  13°. 
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Figure  3. 
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Variations  of  azimuth  wavenumber  location  for  the  half-maximum 
Q2- contour .  Parameter  values  as  given  in  Equation  (54). 
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Figure  4.  Variations  of  range  wavenumber  location  for  the  half -maximum 
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Figure  6.  Surface  wave  propagation  direction  (as  measured  from  the  range 
coordinate)  versus  amplitude  s0  of  the  corrugating  wave. 

The  peak  values  of  Qz  as  given  in  Figure  5  vary  only  over  the  range 
3.35  to  4.90  for  ?o-0.5  metres,  only  1.6  dB,  and  so  the  direct  variation  of 
Q2  is  a  relatively  minor  amount.  For  the  Kelvin  wake  problem,  the  periodic 
variation  of  the  filter  window  in  ky-space  is  potentially  much  more 
significant . 

For  experimental  comparisons,  there  are  three  images  of  clear 
narrow-V  wakes  that  were  taken  during  the  JOWIP  measurements,  those 
designated  JOWIP  8/2,  8/3  and  9/2.  They  are  shown  in  Figure  7.  From  these 
images  the  orientation  of  the  most  steeply  inclined  narrow-V  arms  for  each 
pass  have  been  measured,  with  respect  to  the  SAR  aircraft  track,  and  the 
values  are  given  in  Table  1.  Ship  track  angles  as  determined  from  the 
positioning  equipment10  are  also  included  for  each  of  these  runs. 
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Table  1.  Ship  Track  and  Steepest  Wake  Arm 
Inclination  to  SAR  Track 


SAR  Pass 

Wake  Arm  Angle 

Ship  Track  Angle 

8/2 

18° 

13° 

8/3 

15° 

11° 

9/2 

11° 

9° 

It  can  be  seen  from  the  Table  that  the  wake  arm  inclinations  are  of 
the  same  order  or  are  larger  than  th^.  processor  acceptance  angles  (at  half¬ 
maximum)  for  ?o=0.5  m.  However,  the  variational  effects  just  described  will 
contribute  strongly  toward  enhanced  (or  decreased)  outputs  in  the  SAR  image 
even  for  these  cases. 


Figure  7.  SAR  images  of  narrow-V  wakes 
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4 .  Internal  Waves 


In  order  to  determine  the  effective  impulse  response  function  of 
the  SAR  system  as  it  pertains  to  propagating  internal  wave  modulations  of 
surface  wave  fields,  a  Fourier  description  for  c  will  be  used  in  Equation 
(21)  with  modulational  features  explicitly  incorporated.  The  surface  wave 
field  will  be  specified  using  a  stochastic  description  rather  than  a 
deterministic  one  as  used  in  the  last  section.  Let  ?  be  described  in  terms 
of  linear  dispersion  theory,  i.e.  frequency  a  being  a  given  function  of 
wavenumber  1c ,  but  with  modulating  transform  terms  A: 


s(x',t')  =  JJJ  ps(ic)A(lc,  af  ;x'  ,  t'  )5  (c7-cry-k?U(x'  ,  t  ’)  )eilc*x'  ' iat'  dCdcr 

-oo 

+  ///  ps*( -K)A*( -Ic.CTy  ;x'  ,  t' )5(a+ay-k^U(x'  ,  t '  )  )ei-k*x' ‘  (56) 

-oo 

Here  ps(lc)  is  a  wide  sense  stationary  transform,  A(lc,oy  ;x' ,  t' )  describes  the 
effect  of  the  internal  wave  modulation  on  the  surface  wave  field,  is  the 
horizontal  surface  current  (assumed  to  be  due  to  the  internal  wave  field 
only),  and  cry  is  the  radian  frequency  of  freely  propagating  surface  waves, 
i.e. 


Of  = 


J  gk+Tk^ 


(57) 


In  this  last  equation,  g  is  the  acceleration  of  gravity  and  T  is  the  surface 
tension  of  the  air-water  interface.  In  Equation  (56),  ?  is  guaranteed  to  be 
real  by  the  form  of  the  right-hand  side,  and  the  surface  wave  propagation 
properties  in  the  presence  of  the  advecting  current  field  0  are  accounted 
for  by  the  5 -functions. 


Equation  (56)  can  be  substituted  into  Equation  (21)  with  the  result 


W(xQ)  -  ff  (ps(£){  ///  A(ic,o/;x' ,t')I(x' .Xo.t'Jeik'x'-io/t'-ik'U't'd^.dt^ 


+  Ps*( -«){///  A*(-ic,a/;x'  ,t')I(x'  ,  x0 ,  t '  )ei^x' +ia/t '  -  ik^U' t '  d^«  dt  -  }]d£  (58) 


and  here,  0'  is  an  abbreviated  notation  for  iJ(x',t'). 


The  image  I(xQ)  is  defined  as 

I(x0)  -  |W(x0) |2 


and  so  the  mean  image  is  given  by 


<I(x0)>  -  <|W(x0)|2>  , 


l.e.  , 


<I(xQ)>  -  //  Ps(k)  f .  ./  A(lc,<7y;x'  ,  t'  )A*(ic,ay;x"  ,  t")I(x'  ,x0,  t'  )I*(x"  ,x0,  t") 


+  IS  A(-^,a/;xH,t")A*(-ic,o/;x'  ,t')I(x'  ,x0,t') 


x  I*(x")x0)t")eik*(x'-x")+io/(t'-t")-ikt(U't'-U"t")dx'dt'dx"dt"diJ  (61) 


where  Ps(k)  is  the  power  spectrum  of  the  ps(k)  process,  defined  by 

<Ps(S)Ps*<£')>  ’  Ps (£)*(£-£') 


(62) 
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In  order  to  proceed  further,  it  is  helpful  to  express 
Fourier  transform  and  I  by  a  Fourier-like  transform: 


A(k>/;x^,t')  -  SJJ  elK?x'-iQt^B(2jQ;Sfa/)d^dn 


Q2<^' ,0' ,2";x0)  -  fff  I(x',x0,t')eiK'«x'-iQ't'-iK""U't'd$/dt- 


With  these  substitutions  Equation  (61)  becomes, 


<I(x0»  -  ff  Ps(Ic)| /// 


+  JT  P.(-«)|JJJ  -ic,a/)Q2(ic-^,  -a/+n,£;x0)d£dfi|2d£ 

-  CO  -CO 


and  at  this  stage  it  is  expedient  to  specialize  to  weak  interactions 
that 

A($,CTy;x'  ,t' )  -  1  +  Ai(1c,£7/;x' ,t' )  ;  Ai«l  . 

Therefore , 

<t> B  “  S(t)S( 0)  +  *Ai(2,n;ic,<7/) 

and  the  squared- term  in  the  first  integrand  in  Equation  (65)  becomes 

lQ2(£.<V.£;x0)  +  JIf  ^Al(^.n'.^'<7/)Q2(^+^-£7/fn.^;x0)d^dn|  2 


A  by  a 


(63) 


(64) 


(65) 

so 


(66a, b) 


(67) 
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Because  of  Equation  (66b)  this  becomes 

|Q2(S,ay,ic;x0)  |2  +  Q2(ic,a/,ic;x0)  JJJ  f  a jH-tt 

-00 

+  Q$(lc,cT/,Ie;x0)  JJJ  ^Ai(K,n;ic,ay)Q2(ic+^,ff/-l-n,ic;x0)dKdn 

-00 

and  the  remaining  terms  are  ignorable,  being  0(Ax2)cf  1.  The  squared  term 
in  the  second  integrand  transforms  in  a  similar  manner.  With  this 
substituted  in  Equation  (65),  the  term  in  |  Q2 1 2  can  be  seen  to  be  simply 
<I(x0)>s,  i.e.  the  image  that  would  have  existed  in  the  absence  of  any 
modulation,  and  so,  the  deviation,  due  only  to  the  modulational  effects,  is 
given  by 


<I(xQ)>  -  <I(xQ)>s  -  JJ  Ps(ic)  {Q2(lc,cr/,ic;x0) 

-00 


CO 

x  JJJ  ^Al(^.n.^.£7/)Q$(k+K.ff/+n.^;x0)d^d^)d^ 

-00 


+  JJ  Ps(-£)  {Q2(£»  -a/,E;x0) 


x  JJJ  ^Al(£.n>  -&,<7/)Q?(k?K, -a/+fl,ic;x0)dKdn)dic 

-00 

+  c.c.  (68) 

The  term  Qz  defined  in  Equation  (64)  differs  only  from  the  corresponding 
term  of  Equation  (33)  by  the  presence  of  the  phase  -  shifting  component 
K"*U't'  due  to  the  internal  wave  current.  However  for  typical  current 
values,  i.e.  1  m/s,  and  even  for  large  processor  integration  times,  i.e.  10 
seconds,  the  product  U't'  is  equivalent  to  only  a  10-metre  shift  in 
position.  This  is  to  be  compared  with  a  pixel  size  of  1  to  10  m  in  range 
(the  range  compression  is  already  included  in  I  in  Equation  (64))  or  a 
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beamwidth  of  order  100  m  in  azimuth.  For  Gaussian  beams,  I  is  a  Gaussian 
shape,  and  so  with  the  K‘'?U't'  term  included  up  to  second  order  in  the  x',t' 
variables  (see  equivalent  analysis  in  Section  3,  Equations  (39)  to  (43)),  Q2 
can  be  determined  explicitly  and  it  also  is  a  Gaussian  function  in  K- space. 
Because  Q2  is  essentially  a  Fourier  transform  of  a  Gaussian  function  with 
-10-m  widths  (after  the  t' -integration) ,  it  displays  widths  in  K-space  of 
order  2/10  rad/m  and  these  are  centered  at  its  Bragg  wavenumber,  £B.  Again, 
using  the  numerical  values  of  Equations  (54a-g),  £B  is  19.7  rad/m  along  the 
range  direction  and  zero  along  azimuth,  and  so,  the  bandwidth  of  Q2  in  k- 
space  (as  used  in  Equation  (68))  is  only  -1%  of  its  central  value,  i.e. 
extremely  narrow. 

From  the  interaction  dynamics,  it  can  be  determined  that  <f>kl  is  a 
broad  function  of  £,  and  from  estimates  of  surface  wave  spectra  it  is  known 
that  Ps(£)  is  broad  and  smooth  over  bandwidths  of  order  1%.  These  functions 
can  then  be  approximated  by  their  values  at  £B  and  removed  from  the  k- 
integration  processes  in  Equation  (68)  with  very  little  error.  For  further 
analysis  the  integrands  are  returned  to  Ax*  and  I  using  Equations  (64)  and 
(63),  and  in  the  process  the  £- integration  is  done  after  the  dispersive 
character  of  the  surface  waves  is  incorporated  explicitly  up  to  linear  terms 
in  £  by  the  equation 


<Jf  “  +  (CgB?(k-kB)  .  (69) 

The  £ -  i n t e g r a t i o n  returns  5-functions  with  arguments 
x' -x"+cg(t'  - 1") -tJ'  t'+ft"t" ,  and,  in  principle  these  allow  one  of  the 
x- integrals  to  be  performed.  The  result  is 
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<I(x0)>  -  <I(x0)>s  -  4w2Ps(icB)  ////  I*(xi",x0,t")I(x' ,x0,t')A]*(icB,CTB;x' ,t) 

-CO 


x  e* i{<7B*cgB^B5  (t'  -t")j (^1«  ix'  )dx'dt'dt" 

+  4fr2Ps*(-icB)  JfJJ  I*(x2",x0,t")I(x' ,x0,t')Ai(-^B,aB;x' ,t') 


x  e 


i{aB-CgB^kB) (t' -t")j (x2" ,x' )dx' dt' dt" 


(70) 


with 


xi,2"  -  x'  +  cgB(c'-t")-lJ(x',t')t'+lJ(xi(2",t)"t"  (71) 

implicitly  defining  x"]^  2  (x^  for  upper  sign  and  x2"  for  lower  sign)  .  In 
Equation  (70),  J  is  the  Jacobian  allowing  the  5-functions  to  be  integrated 
and  it  has  the  form 


J 


(1 


t"£Hx)(l 

3x" 


t"£Ex) 

ay" 


evaluated  at  x^  or  x2"  as  appropriate. 


(72) 


The  modulation  Ax  is  usually  described  in  terms  of  the  energy 
envelope  E(£,ay;x' , t' )  in  which  case,  by  defining  the  autocorrelation 
function  of  the  s -process  in  terms  of  the  Fourier  transform  of  Ps(jc)E  and 
equating  that  description  to  the  equivalent  expression  from  Equation  (58), 
the  following  result  pertains: 
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|AOc,a/;x'  ,t')  |2  -  {E+E*}/2 

or,  in  terms  of  At  and  an  equivalent  Et  defined  as  E-l, 

I A  | 2  -  1  +  Ai+Ai*  +  .  . 

-  1  +  (E1+E1*)/2 


and  so , 


(73) 


(74) 


Aj.  “  Ej/2  (75) 

For  usual  applications,  current  strain  rates11  are  less  than 
0(10  2s_1)  and  SAR  processing  times2  are  less  than  10  s,  thus  in  Equations 
(72)  and  (70),  J  can  be  approximated  quite  well  simply  as  unity.  This  also 
simplifies  Equation  (71)  to  the  explicit  form 


xi, 2"  -  x'  +  cgB(t'-t")-U(x' ,t')(t'-f)  (76) 

The  propagational  character  of  the  internal  waves  can  also  be 
incorporated  explicitly,  and  here  for  simplicity,  only  dispersionless 
internal  waves  will  be  considered.  A  single  transformation  can  be  made  from 
x',t'  and  x",t"  frames  to  coordinates  moving  with  the  internal  wave  rj ,  t'  and 
>71,2.  t"  : 


x'  -  rj+cjt' 

(77a) 

xi,  2"  ~  n,2+cit" 

(77b) 

where 

c j  is  the 

phase  (and  group) 

velocity  of  the  internal  waves. 

With 

these 

changes  , 

Equation  (70)  can 

be  rewritten  as  follows  (note, 

from 

Equation  (62),  Ps  is  real): 


-29- 


<I(x0)>  -  <I(x0)>s  -  2w2Ps(icB>  //  Ei*(icg  ,af,rj)  JJ  I(»?+c;[t'  ,x0 ,  t'  ) 


x  I*(^It"^0lt")e-i{£7B-cSB-kB)(t:  -t  )dt'dt”d^ 

+  2jt2Ps(-^b)  //  E]_(  -Sb , <7/ ,  f?)  /J  I(»7+cxt' ,x0t' ) 

-CO  -00 

x  I*^2+SIf,x0,t»)ei{£7B-CgB"kB}(t'-t'-)dt;(dt..d- 

+  complex  conjugate  (78) 

With  a  further  definition  of  the  inner  double  integrals  in  Equation 
(78)  as  impulse  response  functions  Ka>b.  this  last  form  becomes 

<I(xQ)>  -  <I(x0)>s  -  2tt2Ps(Kb)  /J  Ei(icB.^/l^)Ka(^,x0)d^ 


+  2ff2Ps ( -Icb)  JJ  Ei*(-SB,(7/,^)Kb(^,x0)d^ 


+  complex  conjugate 


where 


Ka.bC’J-Xo)  “  J/  I(»?+cit'  ,x0lt')I*(ni,2+cit",x0,t") 


x  e+^B-tgB-kBXt't")^,^,. 


with  Ka  taking  the  upper  sign  and  rji,  and  the  lower  sign  and  rj 2. 
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From  Equation  (79)  it  can  be  seen  that  the  modulation  in  the  image 
is  simply  a  scaled,  filtered  version  of  the  Bragg  scatterer  energy  density 
added  to  a  similar  scaled  filtered  version  of  the  anti-Bragg  energy  density. 
The  kernel  functions  Ka  b  represent  filters  that  depend  on  position  xQ  (and 
are  not  just  functions  of  rj-x0)  . 

For  the  Gaussian  beam  case,  I  is  given  by  Equation  (24),  and,  using 
Equations  (25)  to  (27)  and  after  expanding  the  exponent  in  I  to  second  order 
in  t'  and  t" ,  the  integrals  in  Equation  (80)  can  be  evaluated.  The  results 
can  be  put  into  a  standard  normal  form  as  follows: 


Ka,b^  >xo) 


zV>TG2YG2 

647t2r08c^a2vi2 


ra  bexp[-41n2{(^is)2  +  (?-Y:Zsa,b) 2  } 
ReS 


Res-, 


-  Ba,b(»lx-xs)(fly-ysa,b)]  (81) 

Here,  the  Resx  y-factors  are  full-width  at  half -maximum,  and  the  various 
terms  to  first  order  in  cI/V1,  U/Vj^  and  ( V-Vl)/V1  s.re 


Resx  - 


_ cRqTg _ 

x0(«'2tg4- 1)1/2 


Jln2 


(82) 


Resy  - 


_ 2cRpJ ln2 _ 

YGWc(1+c2Ro2/Yg4wc2)1/2 


X  {1  -  +  V.-Yl  +  j(?.Iy'c.RBy-uy)  ^cYG2 ; 

2Vl  2Vi  V!  cRq 


(83) 


XS  -  Xo-yoCixAl 


(84) 


ysa  -  ySb  -  y0-^y0(v'^-cly)  +  *0  — - 

V1  Vi{1+c2r02/ygV2> 


(85) 
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Ba  -  Bb  -  cgBx+Ux-cIx(2iYG2g,cx02(i-<',2T  4) 
V1  R03TG2c3 


+  £&BxlUx  iflGjzp2  cix  _Jxq2_  (1.  i  2Tg4)  }  (86) 

Vl  c2r04  V!  r02Tg2c2  G 


^cgBv+cIv+^v'i  +  i-ucYc^  ^CgBvcIv+^y^  .  V~Vl 


Fa  **  Fb  ”  e 


2Vl 


cR, 


2Vi 


Vl 


(87) 


The  differences  between  these  equations  and  those  of  the  previous  section 
(see  particularly  Equations  (75)  and  (52)  are  due  completely  to  the  product 
form  of  Equation  (73)  and  its  convolutional  nature.  Resx  here  is  less  than 
Resx  in  the  Kelvin  wake  description  by  a  factor  of  about  $2. 


.  It  should  be  especially  noted  that  the  y -behaviour  of  Ka  b  is  not 

confined  to  the  explicit  fjx,  rjy  terms  appearing  on  the  right-hand  side  of 
Equation  (81).  The  terms  Resy,  ysab  and  Ta  b  also  contain  ^-behaviour 
through  Ux  and  Uy.  The  traditional  velocity  bunching  term  x0Ux/Vt  appears  in 
Equation  (85)  and  it  produces  not  only  the  usual  misregistrations  but  it  can 
also  produce  variations  in  image  intensity  if  there  are  variations  in  Ux 
with  rjy.  This  can  be  seen  by  considering  Ux  to  be  linear  in  r/y  in  the 
exponent  of  Equation  (74)  thus  effectively  changing  the  y-resolution.  In 
symbolic  form,  and  with  c^2 /Yg^w,.2  ignored  compared  to  1  (as  is  typical)  , 
and  U'X0-(3Ux/3y)0, 


ux  “  uxo+»?yUxo 


(88) 


Ys 


yso  + 


Xq>?vUxo 

Vl 


(89) 


V  ,  Z^Wc4TG2yG2  r  -41n?l't?X-  XS)2-£1r.9(tW^-X0U'x0/Vl]  -yS0)2 

a>b  ^r2Ro8c4A2Vl2re  ReS*  Resy 


x  e-B(f?x-xs)('7y[l-^Ya:]-yso) 


(90) 
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If  Ei  were  independent  of  rj ,  the  E^  terms  could  be  removed  from  the 
integrals  in  Equation  (79)  and  the  integrations  performed  using  Equation 
(90).  The  result  would  be 


<I(x0)>  -  <I(x0)>s  -  27r2[Ps(iJB)E1(lcB)+Ps(-icB)E1*(-icB)] 


x  z^Tg2Yg2 

64ff2RQ8c4^2v^2 


?rResxResA 


Resx2Resv2B2 


1 

1-  xqU' o 
VX 


(91) 


with  the  scene  intensity  showing  modulations  due  to  the  last  term  on  the 
right  even  though  no  mean  modulations  were  present  in  the  original  wave 
field.  The  same  analysis  could  have  been  performed  on  the  <I(xD)>s  term 
showing  that  it  was  perturbed  by  the  same  right-hand  term,  and  then  Ex  could 
have  been  set  to  zero  for  the  uniform  wave  field  condition.  Either  way  the 
result  would  be  that  the  gradient  of  Ux  in  y  produces  pseudo -modulations  in 
standard  SAR- processed  imagery. 


For  the  symmetric  processor  case,  yo-0  and  it  can  be  seen  from 
Equations  (82)  to  (87)  that  the  scene  motion  and  aircraft  mismatch  affect 
only  Resy,  ys  and  T,  but  Resx  and  xs  are  unaffected  to  this  approximation. 
For  typical  parameters  (e.g.  Equations  (54a)  to  (54m)) 


“cV  »  i 

cR0 


(92) 


w'Tq2  »  1 


(93) 


and  so 


Resx  -  cR9  Jln2 
xoW'TG 


(94) 


-33- 


Re  sy 


_  a>c2YG4^cIv-cgBv-Uy)2 

2cR0J ln2  1  +  c2rq2  yt _ 

YGWC  X  .  Cly  +  V-Vj  _  4^CYG2  ^■CXy-CgBy-Uy-' 
2V!  2V!  cR0  VL 


(95) 


B  -  ceBx+ux  ^>2Yg2zo2  .  £lx  2x0V2Tg2 
V1  c2Rq^  R02c2 


- 1. 


:gBx 


+ux~cIx 


2YG2q,cx02tj,'2TG2 


Vl 


Ro 


3C3 


(96) 


Xc 


(97) 


Ys  “  x0.cRBxt-Hx  (9g) 

and  T  is  unchanged.  With  Resy  written  as  in  Equation  (95),  the  exponent  of 
Ka,b  from  Equation  (81)  can  be  seen  to  separate  quite  readily  into  a  real 
and  an  imaginary  part.  The  resultant  real  resolution  in  y  is  thus  obtained 
from  the  exponential  of  the  real  part  of  the  Ka  b  exponent,  multiplied  by 
the  cosine  of  the  imaginary  part.  The  resolution  width  occurs  when  this 
product  equals  1/2.  In  symbolic  form,  and  with  the  Ka  b  exponent  defined  to 
be  the  resolution  width  is  obtained  from 

e^rcos«^x  -  1/2  .  (99) 

By  using  Equations  (90)  and  (95)  at  r;x-xs  and  with  some  algebraic 
manipulations,  Equation  (99)  provides  the  result 


(Resy)Real 


“t0*  ^i(i  -  a? +  ^>2  +  <«»> 


-  (i 


cIv 

2Vl 


)2  - 


r%  / 


where 


a 


«czYG4(CIv-CgBv-Uv)2 

c2rq2  Vx 


(101) 


a/b  ■+  -hd 
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and 

(Resy) real  -  2 ^  cR°  .  (105d) 

This  is  the  traditional  value  and  it  decreases  as  the  width  of  the  processed 
beam,  YG,  increases. 

For  cIy,  cgB,  U,  *0  and  Yg-"*j, 

a/b  -*  -oo  (106a) 

9  -  k/ 2  (106b) 

and 


(Resy)real  -  YG2E  |  cIy'cgBy-%  1 3/2  .  (106c) 

CK0  Vi 

In  this  limit  the  resolution  increases  as  the  beamwidth  YG  increases  and  it 
depends  on  the  azimuthal  translation  speed  of  the  surface  wave  Bragg 
scatterers  (cgBy+Uy)  relative  to  the  azimuthal  component  of  the  internal  wave 
phase  velocity  (cIy)  . 

For  the  numerical  values  of  Equations  (54a)  to  (54h) ,  and  with 

ciy  -  (CgBy+Uy)  -  1  m/s  (107) 

the  results  from  Equations  (99)  to  (101)  are 


a/b  -  0.4661 

(108a) 

9  -  0.785 

(108b) 

(Resy)real  =>  4.19  m 

(109) 
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This  can  be  compared  with  3.16  m  obtained  from  Equation  (105d)  showing  that 
a  not  unreasonable  amount  of  azimuthal  scene  motion  can  increase  the 
resolution  size  by  a  significant  amount  (33%  in  this  case). 

To  complete  the  description  of  the  effect  of  scene  motion  it  can  be 
seen  from  Equation  (98)  that  the  velocity  bunching  term  should  actually  be 
the  range -component  of  the  translational  velocity  of  the  Bragg-scatterers 
cgBx  +  rather  than  just  Ux  itself.  Furthermore  there  are  small  amplitude 
and  phase  variations  of  the  K's  due  to  the  azimuthal  components  (and  the 
processor  mismatch)  as  is  evident  in  the  T's  (Equation  (87)).  The  amplitude 
terms  are  only  of  order  1-10%  and  the  phase  term  is  of  order  1/2  radian. 
There  is  also  a  small  eccentricity  of  the  impulse  response  peak  as  shown  by 
the  non-zero  value  of  Ba,b  in  Equation  (86),  and  for  typical  values  this  is 
of  order  0.40  and  dependent  solely  on  the  range  -  components  of  the  scene 
velocities . 

From  Equation  (59)  the  covariance  of  the  image  can  be  determined, 
and  using  the  assumption  of  Gaussian  statistics  for  the  (narrow-band) 
transform  coefficients  ps(&)  ,  the  covariance  can  be  shown  to  be  (see 
Appendix) 


covj  (A)  -  16jt4Ps2(Icb)  [  | //{ Ipl (x , x0) I$i(x , x0+2)  +  Ip2(x,x0)I$2(x,x0+A)  }dx|  2 


00 

x  | J/{Ipi(x,x0)Ip2(x,x0+2)  +  I pi (x , x0+A) Ip2 (x, xD) ) dx | 2 ]  (110) 

-CO 

where  IP1  is  the  unsquared  impulse  response  function  (i.e.  whose  width 
defines  the  pixel  size  of  w(x0))  for  positive  of  and  IP2  is  the  same  for 
negative  of.  It  can  be  seen  that  if  \~&\  is  much  larger  than  one  pixel 
dimension  the  covariance  is  essentially  zero  because  the  products  of  the 
IP's  are  then  zero. 
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5 .  Discussion 


Kelvin  wake  narrow-V's  in  the  SAR  imagery  are  considered  to  be 
partly  due  to  some  randomization  and  intermittency  in  the  wake  structure 
itself,  as  is  inferred  from  experimental  evidence12,  however,  the 
peculiarities  of  the  SAR  processor  described  above  do  provide  an  alternative 
mechanism  for  enhancement  of  local  regions  of  the  wake.  The  fact  that  these 
are  predicted  to  occur  once  per  Cycle  of  the  transverse  wave,  and  for  the 
enhancements  in  each  of  the  two  arms  of  the  narrow-V  to  occur  at  the  same 
transverse  wave  phase  angle,  also  lends  credence  to  the  validity  of  the 
effect,  because  these  are  in  accordance  with  experimental  data.  No  attempt 
has  been  made  to  model  fully  the  Kelvin  wake  field,  transverse  and 
diverging,  because  it  is  known  that  accurate  results  for  the  short 
wavelength  components  are  very  difficult  to  achieve,  requiring  complicated 
source  descriptions  at  the  vicinity  of  the  source  ship,  and  also  because  of 
the  intermittency  just  referred  to.  Indeed,  the  predicted  SAR  processor 
variations  are  not  large  enough  by  themselves,  but  combined  with  the 
intermittency  and  randomization  of  the  direction  of  the  diverging  wavefield, 
they  could  possibly  provide  the  mechanism  underlying  the  appearance  of  the 
brightened  V-arms. 

The  curvature  effect  examined  in  Section  3,  Equations  (50)  to  (53), 
can  be  much  larger  than  traditional  velocity  bunching  effects.  The  latter 
are  known  to  produce  radar  cross  section  aberrations  that  are  given  by9 

KbI  -  za\?hi\  sa  dll) 

dyz  V1 

where  the  pattern  is  considered  to  be  moving  in  the  y-direction  as  a  solid 
body.  The  comparable  term  from  the  curvature  effect  is 

kcurl  -  2<SOS£s)2  (112) 

dy 2  wc^G 
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which  is  obtained  by  integrating  K(rj,x0)  over  all  rj ,  as  if  the  Bragg - 
scatterers  covered  the  entire  pixel  area  uniformly,  and  using  Equations 
(51),  (52)  and  (53).  In  Equations  (111)  and  (112)  only  the  perturbations  in 
a  are  given.  The  denominator  in  Resy  from  Equation  (53)  is  also  expanded 
binomially  and  the  first  non-unity  term  kept. 

For  the  numerical  values  of  Equations  (54a-j),  the  relative 
magnitudes  of  cryj  and  dcur  are  0.06  and  0.5,  each  of  these  values  coming  from 
the  last  terms  in  Equations  (111)  and  (112)  respectively.  This  evaluation 
is  only  for  a  very  specific  kind  of  surface  (parabolic  within  the  pixel 
area)  and  so  is  only  indicative  of  possibilities.  But  it  does  show  that 
curvature  effects  are  not  ignorable,  a  priori. 

The  perturbation  in  the  SAR  processor  output  due  to  the  curvature 
of  the  corrugating  wave  can  be  understood  from  a  simple  geometrical  picture 
of  the  radar  wave  at  the  sea  surface.  In  the  absence  of  the  corrugation,  a 
field  of  plane  radar  waves  incident  on  the  sea  surface  at  an  angle  0inc 
intersect  the  surface  in  a  series  of  straight  lines  as  shown  in  Figure  9(a). 
In  the  presence  of  a  corrugation  propagating  perpendicularly  to  the  radar 
waves  (in  the  horizontal),  the  intersection  of  the  radar  waves  at  the 
surface  also  becomes  corrugated,  with  local  horizontal  wavenumber  directions 
that  can  be  considerably  deflected  azimuthally  from  that  of  the  main  radar 
beam,  Figure  9(b). 

From  Figure  10(a),  it  can  be  seen  that,  with  z  upwards,  the  radar 
crest  that  just  intersects  the  sea  surface,  as  shown,  does  so  at 

zt  -  (xr-ztanfltne) tan0inc  (113) 
If  the  surface  height  is  specified  by  a  long  surface  wave  with  height  ? 
given  by 

?  -  ?Qcos  (k?x+<£)  (114) 
where  s0  is  the  wave  amplitude,  Ic  is  its  wavenumber,  x-(xr,xa)  are  range  and 
azimuth  coordinates,  and  4>  is  an  arbitrary  phase  angle,  then  the 
intersection  takes  place  at 


Zi  -  T 


(115) 


which,  by  Equations  (113)  and  (114)  is  at 


xr-ztan0inc  -  — is — cos(k?x+<£)  .  (116) 

tannine 


Surface 


Figure  9.  Intersection  of  the  radar  beam  and  the  ocean  surface  (a)  non- 
corrugated,  (b)  corrugated. 

This  last  can  be  solved  for  xt  in  terms  of  x,  (in  principle)  giving  the 
horizontal  shape  of  the  intersection  lines.  For  a  corrugation  propagating 
along  the  azimuth,  k?x  -  k,x,  only,  and 


xr  -  — is - cos (kaxa - <f> )  +z tanS inc  .  (117) 

tanS inc 

Typical  plan  view  intersection  lines  are  shown  for  this  case  in  Figure 
10(b). 
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Figure  10.  (a)  Geometry  of  intersection  (elevation  view).  (b)Typical 
horizontal  intersection  lines  for  a  radar  field  in  the  presence 
of  an  azimuthally  propagating  corrugation  (plan  view) . 

The  slopes  of  these  lines  provide  estimates  of  the  azimuthal 
deflection  amounts  for  the  effective  radar  beam  in  the  Bragg  scattering 
process.  These  deflections  allow  the  scattering  to  take  place  from  a  wider 
range  of  horizontal  surface  wave  propagation  directions  (in  the  azimuth) 
than  would  otherwise  be  possible.  From  Equation  (117)  the  slopes  are 


2*?q 

Astan0 inc 


+  <f>) 


(118) 


where  As  is  the  wavelength.  For  the  values  of  Equation  (54a-m)  but  with 
?o  -  0.5  m, 

|^Sl|  -  0.177  -  tan  10° 
dxa 


(119) 
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This  is  in  accordance  with  the  results  of  Section  3  where,  from  Figure  6, 
the  centre  of  the  SAR  processor  at  0-0°  is  aligned  with  waves  propagating  at 
9.6°  (at  the  maximum  deviation). 

The  internal  wave  problem  has  been  modelled  by  Fourier  transforms 
and  linear  dispersion  theory  for  the  surface  waves.  This  essentially 
resolves  the  wavefield  into  sinusoids  of  uniform  height  but  of  random 
amplitudes  and  phases  (across  different  realizations) .  The  measure  of  the 
coherence  of  the  wave  field  is  contained  in  its  power  spectrum  and  the 
bandwidths  imposed  by  the  radar  configuration  and  the  processor.  In  all 
cases,  the  surface  waves  imaged  by  the  SAR/Bragg- scatter  process  are 
coherent  over  essentially  one  pixel  dimension,  regardless  of  the  pixel  size 
(for  pixels  large  compared  to  a  Bragg  wavelength) .  This  is  apparent  for  the 
present  model  from  Equation  (110)  and  it  follows  directly  from  the  fact  that 
for  larger  pixels  i.e.  wider  behaviour  in  Ip,  better  resolution  is  obtained 
in  the  Fourier  domain,  i.e.  smaller  bandwidths  pertain  and  Q2  is  narrower, 
but  these  result  directly  in  larger  coherence  distances.  In  fact  the  two 
are  quite  simply  complementary  though  the  narrow-bandedness  of  the  Fourier 
description  of  the  power  spectrum.  For  this  particular  model,  broadband 
coherence  distances  are  irrelevant;  only  the  narrow-band  coherences  enter 
the  problem.  However,  if  strong  non-linearities  were  present  in  the  surface 
waves  then  this  simple  conclusion  would  need  to  be  reexamined. 

At  each  pixel  location  the  mean  image  (Equation  (60)  rather  than 
Equation  (59)),  because  of  its  squared  nature  and  all  of  the  SAR  processing, 
is  due  ultimately  to  the  power  in  the  Bragg  waves  covering  the  entire  pixel 
area.  If  the  Bragg  wavelengths  were  much  smaller  than  a  pixel  dimension  so 
that  individual  scattering  centers  could  be  postulated  within  the  pixel 
area,  then  internal  correlations  between  Bragg  waves  resolved  around  their 
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individual  centres  could  be  established,  and,  in  keeping  with  previous 
work,1  would  likely  show  uncorrelatedness .  In  that  case  the  bandwidth  of 
the  scattering  process  would  need  to  be  large  enough  to  allow  that  detailed 
structure  to  be  faithfully  captured  in  the  scattered  radar  beam,  and  the  SAR 
processor  would  need  to  be  wide  enough  in  bandwidth  to  accept  it,  otherwise 
the  internal  structure  within  the  pixel  would  be  smoothed  out  and  some  of 
the  uncorrelatedness  would  disappear.  When  the  bandwidth  of  the  processor 
is  due  to  the  pixel  dimensions,  as  is  the  case  here  (Equations  (64)  or 
A(6)),  essentially  no  structure  is  allowed  in  the  scattering  waves  within 
the  pixel,  and  so  internal  decorrelation  distances  (or  times)  cease  to  be 
meaningful,  and  the  signal  reflected  from  the  pixel  becomes  uniform  across 
the  pixel  area,  i.e.  internally  coherent.  The  power  in  this  coherent 
signal,  as  stated  earlier,  is  equal  to  the  average  power  over  the  pixel 
area.  This  occurs  because  the  Bragg- scattering  process  is  essentially  a 
Fourier  transforming  process  with  the  SAR  processor  setting  the  bandwidths 
and  these  in  turn  setting  the  pixel  size.  It  should  be  noted  that  this 
concept  has  been  explored  previously  in  connection  with  surface  waves.2 

The  variation  of  Resy  with  YG  (Equation  (99))  is  highly  significant 
in  attempts  to  produce  very  small  pixel  SAR's  to  image  internal  waves.  In 
the  example  shown  in  the  last  section,  Equations  (106)  to  (109),  a  processor 
that  produces  a  Resy-value  of  3.16  m  in  the  absence  of  azimuth  velocities 
provides  a  value  of  4.19  m  in  the  presence  of  1  m/s  components  overall  (see 
Equation  (107)).  If  YG  is  increased  to  256  m  so  that  the  unperturbed  pixel 
size  decreases  to  2.0  m,  the  effect  of  the  azimuthal  components  is  to 
increase  this  resolution  width  to  5.4  m,  and  in  fact,  further  increases  in 
Yg  result  in  even  larger  resolution  widths.  For  these  parameters  an  L-band 
Resy-value  of  2  m  is  possible  but  only  if  R0  is  560  km  for  YG  at  162.2  m, 
and  an  even  larger  Rq  for  YG-256  m. 

For  currents  and  internal  waves  moving  in  the  x-direction  only, 
Resx  and  Resy  both  take  on  their  unperturbed  values  but  the  processor  peak 
is  rotated  by  B,  and  shifted  by  y,  (Equations  (96),  (98)).  The  B-value  also 
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produces  some  eccentricity,  and  if  the  unperturbed  processor  Resx-  and  Resy- 
values  are  equal  (as  they  are  for  the  parameters  of  Section  3)  ,  the 
eccentricity  is  of  order  (Ux/V  i)1/2 ,  i.e.  -0.1  and  so  not  of  great 
importance.  Here,  Ux  is  the  x-component  of  the  surface  velocity,  and 
eccentricity  is  defined  as  for  an  ellipse,  i.e.  J  a2  -b2/a  where  a  and  b  are 
the  semi -major  axes  respectively. 

6 .  Conclusion 


For  the  Kelvin  wake  problem  the  main  conclusion  is  that  large  scale 
surface  curvatures  can  exert  a  major  influence  on  the  processed  SAR  imagery. 
For  particular  cases,  as  examined  in  Section  3,  the  curvatures  associated 
with  the  corrugating  effect  of  the  wake's  transverse  wave  component  can 
readily  rotate  the  processor's  Bragg  acceptance  angle  from  being  bore- 
sighted  with  the  radar  beam  to  being  -10°  fore  or  aft  of  the  beam.  In 
modelling  SAR  imaging  of  Kelvin  wakes,  particularly  for  narrow-V  effects, 
this  mechanism  should  be  included. 

For  internal  waves,  the  main  conclusion  is  that  azimuthal 
propagational  velocities  can  strongly  affect  the  resultant  azimuthal 
resolution  width  of  the  processor.  Also,  because  of  the  Fourier  description 
and  the  use  of  linear  dispersion  theory  for  surface  waves,  and  lst-order 
Bragg  scattering  for  the  radar  reflection,  the  decorrelation  distances  are 
all  of  order  one  pixel  size.  Velocity  bunching  is  shown  to  be  caused  by 
CgBx+Ujc  rather  than  just  Ux,  where  CgBx  is  the  range  component  of  the  Bragg 
wave  group  velocity  and  Ux  is  the  range -component  of  the  horizontal  surface 
current  (vertical  velocities  have  been  specifically  excluded  for  this 
problem) .  Velocity  bunching  also  produces  pseudo-modulation  in  the  imaging 
if  3Ux/3y  is  non- zero,  where  y  is  the  azimuth  coordinate. 
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Appendix  A 


Image  Covariance 

Equation  (58)  can  be  written  as 


w(xQ)  “  JT  ps(£)Kl(ic,x0)dic  +  //  p*(ic)K2(ic,x0)dic  (Al) 

-CO  -00 

where  Kx  and  K2  represent  the  triple  integrals  of  the  A-  and  I -functions. 
The  covariance  of  the  image  is  defined  as 

covj(A)  -  <I(x0+A)I(x0)>  -  <I(x0+£)>  <I(x0)>  (A2) 

-  <|w(x0+A) l2|w(x0) |2>  -  <|w(x0+A) |2>  <|w(xQ) |2>  (A3) 

The  latter  form  of  this  equation  can  be  written  out  in  full  using  Equation 
(41) ,  and  the  expectation  operator  when  applied  to  the  quadruple  products  of 
the  ps(ic)'s  can  be  seen  to  return  double  products  of  the  power  spectrum 
P3(lc),  or  zero  (using  Gaussian  statistics  for  ps(£)).  After  some  lengthy 
algebraic  analysis  using  Equations  (63),  (64)  and  (67),  the  result  for  covj 
becomes 


covj  (A)  -  |/J  Ps(lc)  (Kj(ic,x0)K*(ic,x0+A)+K2(ic,x0)K2(^,x0+^)  }d£|2 

-  oo 


+  |J7  Ps(ic){Kj(ic,x0)K2(ic,Xo+A)+Kj(ic,x0+A)K2(^,x0))d^|2,  (A4) 

-00 


and  on  keeping  only  the  largest  terms  Equation  (A4)  becomes, 


-A2- 


0°  ^ 
covj(^)  -  |JJ  Pg(S)  {Q2(ic,<7/,ic;x0)Q5(ic,CT/,lc;x0+2) 
-  00 


+  Q2(£,  -£7y,Ic;x0)Q5(^,  -ct/,1 c;xq+^)  > die |  2 

+  I//  Ps(£)  {Q2<^.<7/.^‘>^o)Q2^. -<7/.^;x0+2) 
-00 


+  Q2(it,ay,S;x0+A)Q2(S, -<7/,ic,x0)  }dlc|  2  (A5) 

By  using  the  argument  given  in  the  main  text  that  Q2  is  a  narrow- 
band  function  of  Ic ,  it  is  possible  to  remove  the  Ps-function  from  the 
integrals.  Also,  the  function  Q2  is  the  Fourier  transform  of  the  spatial 
impulse  response  function  of  w(xQ)  ,  i.e.  the  smallest  resolvable  element  in 
the  imaged  scene  (unsquared) ,  and  this  can  be  incorporated  in  the  previous 
equation  using  the  following  definition 


Q2(ic,±<7/,£,Xo)  -  JJ  eik*xIplj2(x,x0)dx  (A6) 

-CO 

where  +of  selects  the  impulse  response  function  Ipl  and  -of  selects  Ip2. 
With  this  definition,  and  the  removal  of  the  Ps(£)'s  from  the  integration, 
Equation  (A5)  becomes  Equation  (110)  of  the  main  text. 

The  impulse  response  nature  of  Ipi>2  can  be  seen  from  Equation  (64) 
and  the  Fourier  transform  inverse  to  Equation  (A6) ,  which  provide 


Ipl,2<x,x0)  -  _L_  /Jj*  I(x',x0,t') 

4ff^  .00 


J*J  e^-k*  (x'  -x)+ic7y(k)  t'  -  ik*U' t '  diJdx'  dt  (A7) 


-A3- 


and,  by  using  Equation  (69)  the  inner  double  integral  can  be  performed, 
resulting  in  a  8 -function.  By  using  the  Jacobian  agrument  leading  up  to 
Equation  (76),  the  5-function  can  be  integrated  and  the  final  result  is 


Ipl,2(x,x0)  -J  I(x±cgBt'+lJ(x)t' ,x0, t' )e+i(aB-cgB*kB)t'dt'  (A8) 

-  00 

These  can  be  identified  directly  from  Equation  (21)  as  impulse  response 
functions  of  w(xQ)  ,  if  the  impulse  propagates  with  velocity  ±cgB+fl  and 
displays  phase  shifts  of  ±(aB-cgB?kB)t'  ,  i.e.  if 


<T(x',t')  -  <S(x' -[x±cgBt'+0(x)t'  ])e+i(crB-cgB*kB)t'  (A9) 

This  is  precisely  in  accordance  with  the  linear  dispersive  nature  of  the 
surface  waves . 
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